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1. Introduction 



Systems brought rapidly out of an initial state into a region in parameter space which 
is characterized by several competing stationary states may undergo ageing behaviour. 
Dynamical scaling and universality was first noticed in the mechanical properties of 
several glass-forming systems rapidly quenched into their glassy phase pQ and has since 
been found and studied intensively in a large variety of systems relaxing towards an 
equilibrium state, see [21 EH El El IE] for recent reviews. In what follows we shall 
restrict to systems without any macroscopic conservation laws and to systems without 
frustrations, conditions which are paradigmatically met in simple ferromagnets with 
dynamics which satisfy detailed balance. 

Convenient tools for the study of ageing behaviour in such systems |0] are the two- 
time autocorrelation and autoresponse functions, which in the ageing regime t,s > 1 
and t — s ^> 1 are expected to show the scaling behaviour 

C(t,s) = (<j ) (t) ( j ) (s)} ~s- b f c (t/s) (1) 

~ s-^Mt/s) (2) 

h=0 

Here <f>(t) is the order parameter and h(s) is the conjugate field, t is called the observation 
time and s the waiting time. For large arguments y — *■ oo, one generically expects 



R(t,s) 



5h(s) 



fc(y) ~ y- Xc/z , fn(y) ~ y~ XR/z (3) 

where Ac and Xr, respectively, are known as autocorrelation [TQ1 HI] and autoresponse 
exponents While the ageing exponents a and b can be expressed in terms of the 
dynamical exponent z and equilibrium exponents, the exponents \c,r are independent 
of these but related to the so-called initial slip exponents [IS]. For critical quenches 



20 

a = b = — , at T = T c (4) 

vz 

(where and v are the usual equilibrium critical exponents) is a consequence of the 
fluctuation-dissipation theorem and of time-translation invariance in the scale-invariant 
equilibrium steady-state. 

Turning to the scaling functions, it has been suggested that their form 

could be determined from the requirement of covariance under so-called local scale- 
transformations which are constructed from the requirement of including the special 
conformal transformations t i— > (at + 0)/ (jt + 5) with a5 — 0^ = 1 in time. For the 
response function this leads to [HI [121 EB] 

R{t,s) = s~ l - a f R {t/s) , f R (y) = foy 1+a '- XR/z (y-i)- 1 - a ' (5) 



where a' is a new independent exponent and fo is a normalization constant. This form 
(or equivalently an integrated response) describes very well all available numerical data 
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for ferromagnetic systems quenched to T < T c [T3J [TSJ HZ] and was shown to be exact 
in several exactly solvable models, see [3] and refs. therein. In these systems, z = 2 
On the other hand, for quenches to T — T c , one has in general z ^ 2, but the agreement 
of (jSJ) with numerical data is still almost perfect [UJ HB1 HH] • However, a second-order 
e-expansion calculation in the 0(n)-symmetric 4 -field theory produced a small but 
systematic correction with respect to [5] and there is support of this finding from a 
numerical study of R(t, s) in momentum space [201 • While in most systems studied so 
far a = a' holds true, several models with a ^ a' are also known [TH] . 

The approach of local scale-invariance uses the dynamical symmetry of 
deterministic equations, whereas the influence of noise is essential in the understanding 
of non-equilibrium dynamics. However, considering a description of the ageing system 
in terms of a Langevin equation, it can be shown that if 'the deterministic (i.e. noiseless) 
part satisfies local scale-invariance, then (i) the response function R(t, s) is noise- 
independent and (ii) the autocorrelator C(t, s) can be reduced to noise- less three- and 
four-point response functions |21j . A further extension of local scale-transformation to 
include the diffusion constant as a new dynamical variable then permits to determine 
C(t, s) and the result was found to be in good agreement with simulational data in the 
2D Ising model quenched to T < T c [2*2*j . 

An important ingredient in the ageing studies discussed so far is the assumption 
of detailed balance for the dynamics. This begs the question what might happen if 
that condition is relaxed. Indeed, numerical studies of the contact process, the simplest 
system of that kind, gave the following results l2~^j: 

(i) Dynamical scaling and ageing only occur at the critical point. This is expected 
since both inside the active and the inactive phases there merely is a single stable 
stationary state. 

(ii) At criticality, the scaling forms eqs. (II 1211 hold true for the (connected) 
autocorrelator and the response function, but with the scaling relation 

a + l = b=^-, (6) 

v±z 

in contrast to eq. with f3 and u± being now standard steady-state exponents. 

In order to get a better understanding of these results, it would be helpful to study the 
ageing behaviour in exactly solvable but non-trivial models without detailed balance. 
Remarkably, it has been realized by Houchmandzadeh and by Paessens and Schiitz 
[2*o] that the bosonic versions of the contact process and of the critical pair-contact 
process (where an arbitrary number of particles are allowed on each lattice site) are 
exactly solvable, at least to the extent that the dynamical scaling behaviour of equal- 
time correlators can be analysed exactly [2HJ 12Z] ■ Here we extend their work by 
means of an exact calculation of the two-time correlation and response functions for 
the bosonic contact process and the critical bosonic pair-contact process. In section 2, 
we define the models and write down the closed systems of equations of motion for the 
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correlation and response functions. We also recall the existing results on the single- 
time correlators I n section 3, we discuss the bosonic contact process and in 
section 4, we describe our results for the critical pair-contact process. In section 5 we 
give the results for the two-time response functions. As we shall see, the critical bosonic 
pair-contact process provides a further example of a model where a and b are different. 
In section 6 we conclude. A detailed discussion of local scale-invariance in these models 
will be presented in a sequel paper. 

2. The models 

Consider the following stochastic process: on an infinite <i-dimensional hypercubic lattice 
particles move diffusively with rate D in each spatial direction. Each site may contain an 
arbitrary non-negative number of particles. Furthermore, on any given site the following 
reactions for the particles A are allowed 



It is to be understood that on a given site, out of any set of m particles k additional 
particles are created with rate \i and £ particles are destroyed out of any set of p > £ 
particles with rate A. Diffusion applies on single particles. We shall be concerned with 
two special cases: 

(i) the bosonic contact process, where p = m = 1, hence £ — 1. The value of k is 
unimportant and will be fixed to k — 1 as well. 

(ii) the bosonic pair-contact process, where p = m = 2. 

While the bosonic contact process arose from a study on the origin of clustering 
in biology [2H], the bosonic pair-contact process as defined here j2B] is an offshoot of 
a continuing debate about the critical behaviour of the diffusive pair-contact process 
(PCPD), see [2E] for a recent review. Initially, this model was introduced [21] in an 
attempt to understand the meaning of 'imaginary' versus 'real' noise but the associated 
field theory turned out to be unrenormalizable fl§\ I3*U] . A lattice version (with the 
'fermionic' constraint of not more than one particle per site) of the model contains the 
reactions 2A — > and 2A — > 3A together with single-particle diffusion A0 <-> 0/1 and 
was first studied numerically in [21]. An intense debate on the universality class of 
this model followed, see [2E], and several mutually exclusive conclusions on the critical 
behaviour continue to be drawn, see [3UJ [321 [SHI EH EH] for recent work. The bosonic 
pair-contact process has a dynamic exponent z = 2 [2E] and is hence distinct from 
the PCPD where z < 2. Its study will not so much shed light on any open question 
concerning the PCPD but it should rather be viewed as a non-trivial example of an 
exactly solvable non-equilibrium many-body system to be studied in its own right. 




pA 



(p — £)A ; with rate A 



(7) 
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The master equation is written in a quantum hamiltonian formulation as d t \P{t)) = 
—H \P(t)) |37| I3H] where \P(t)) is the time- dependent state vector and the hamiltonian 
H can be expressed in terms of annihilation and creation operators a(x) and af(x). We 
define also the particle number operator as n(x) = a\x)a{x). Then the Hamiltonian 
of the model (JZJ) reads j2Sj 

d 

H = — D \a(x)a)(x + e r ) + gt(a;)a(a; + e r ) — 2n(a;) 



r=l a; 



- A 



X 



(a\x)) P e (a(x)) p - Y[(n(x)-i + l) 



i=i 



(8) 



(at(a;)) m+/c (afc))" 1 - JJ (n(aj) -2 + 1) 



i=i 



where e r is the r th unit vector. For later use in the calculation of response functions we 
have also added an external field which describes the spontaneous creation of a single 
particle — > A with a site-dependent rate h = h(x, t) on the site x. 

Single-time observables g(x,t) can be obtained from the time-independent 
quantities g(x) by switching to the Heisenberg picture. They satisfy the usual 
Heisenberg equation of motion, from which the differential equations for the desired 
quantities can be obtained. The space-time-dependent particle-density p(x,t) : = 
(a\x,t)a(x,t)) = (a(x,t)) satisfies 



0_ 

dt 



(a(x, t)) = DA X (a(x, t)) - \t (a(x, t) p ) + fik (a(x, t) m ) + h(x, t) (9) 



where we have used the short-hand 

d 



A x f(x) := (/(* " ^, t) + f{x + e r , t) - 2f(x, t)) 



(10) 



=i 



and similar equations hold for the equal-time two-point correlation functions, see 
It turns out that for the bosonic contact process p = m = 1 these equations close for 
arbitrary values of the rates. On the other hand, for the bosonic pair-contact process 
where p = m = 2a closed system of equations is only found along the critical line given 
by 



i\ = pk. 



This line separates an active phase with a formally infinite particle-density in the 
steady-state from an absorbing phase where the steady-state particle-density vanishes, 
see figure Q for the schematic phase- diagrams. In what follows, the essential control 
parameter is 



a := /ik(k + £)/(2D) 



(12) 
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Figure 1. Schematic phase-diagrams for J) / of (a) the bosonic contact process 
and the bosonic pair-contact process in d < 2 dimensions and (b) the bosonic pair- 
contact process in d > 2 dimensions. The active region 1, where lim^oo p(x, t) = 0, is 
separated by the critical line eq. from the absorbing region 2, where p(x,t) — > oo 
as t — > oo. On the critical line p(i) := J da; p(a;, t) remains constant. By varying a one 
moves along the critical line. Along the critical line, one may have clustering (full lines 
in (a) and (b)), but in the bosonic pair-contact process with d > 2 the steady-state 
may also be homogeneous (broken line in (b)). These two regimes are separated by a 
multicritical point. 



The physical nature of this transition becomes apparent when equal-time correlations are 
studied [23 and can be formulated in terms of a clustering transition. By clustering 
we mean that particles accumulate on very few lattice sites while the other ones remain 
empty. Now, for the bosonic contact process, the behaviour along the critical line is 
independent of a. If d < 2, there is always clustering, while there is no clustering for 
d > 2. On the other hand, in the bosonic pair-contact process, there is on the critical 
line a multicritical point at a = ac, with 

1 f°° 

a c = ac(d) = — , A 1 := du {e~ Au I (Au)) d (13) 
^i Jo 

and where Io(u) is a modified Bessel function such that clustering occurs for a > a c 
only and with a more or less homogeneous state for a < ac- Specific values are 
«c*(3) ~ 3.99 and ac(4) ~ 6.45 and limrf\ 2 «c(^) = 0. We are interested in studying 
the impact of this clustering transition on the two-time correlations and linear responses. 

In order to obtain the equations of motion of the two-time correlator, the time- 
ordering of the operators a(x, t) must be taken in account. From the Hamiltonian 
eq. (jHJ) without an external field h, we get the following equations of motion for the 
two-time correlator, after rescaling the times t i— > t/(2D), s i— > s/(2D), and for t > s, 
|4*0] (for a detailed computation, see |4*T| ) 
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d 

— (a(x,t)a(y,s)) (14) 

1 \i pk 

= -A x (a(x, t)a(y, s)) - — (a(x, t) p a(y, s)) + — (a(x, t) m a(y, s)) 

which we are going to study in the next sections. 
3. The bosonic contact process 

For the bosonic contact process, we have p = m = 1, hence also £ = k = 1. We first 
consider the critical case X£ = pk. We shall assume throughout that spatial translation- 
invariance holds and use the notation 

F(r; t, s) := (a(x, i)a(x + r, s)) . (15) 

Then F satisfies a diffusion equation which is solved in a standard way by Fourier 
transforms. It is easy to see that the solution of the equations of motion (|14|) involves 
the single-time correlator F(r,t) := F(r;t,t) which satisfies the equation of motion, 
after the usual rescaling t i— > t/(2D) eq. (10)], 

.9 

— F(r, t) = A r F(r, t) + ap 5 r , (16) 

and the parameter a was defined in (112)1 . As initial conditions, we shall use throughout 
the Poisson distribution F(r, 0) = p\. Hence one arrives at the following expression of 
our main quantity of interest, the connected correlator f 

G(r; t, s) := F(r; t, s)-p 2 = ap dr b (r, ^(t + s) - rj (17) 

where {I r {t) being a modified Bessel function) 

b(r,t) =e- 2dt I ri (2t)...I rd (2t). (18) 

We evaluate this expression in two cases 

• r — 0, t and s in the ageing regime: 

In this case both s and t — s are large, so that we can use the asymptotic behaviour 
I (t) ~ (27rt)" 1//2 e t for t large |12] for the expression 6(0, |(i + s) — r) under the 
integral in (|T7)l . We have to distinguish the cases d > 2 , ol = 2 and d < 2. For 
d > 2 we obtain 

f In 12] this same quantity was denoted by T(t, s) which we avoid here in order not to create 
confusion with the incomplete gamma function 39 . 
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(47T)i(f-l 



-#+1 fj- , „\ -f+1' 



(19) 



By analogy with eqs. (JTJ) and (JHJ), we expect the scaling behaviour G(t,s) : = 
G(0; t, s) = s~ b fa{t/s). We read off the value b = ~ — 1 and the scaling function 

= ip^I) (C - I)"** 1 - to + ■ <*» 

From the expected asymptotics /gO/) ~ y~ x °/ z for y ^> 1, we obtain 

Ac = d (21) 

as can be seen from the asymptotic development of ({20]) and where we anticipated 
that the dynamical exponent z = 2, see also [22] and below. 

For al = 2 the integral in (J 19)) gives a different result. We find 

G(t, s) = f G (t/s) , f G (y) = In (*±T) (22) 

2(2vr)2 \2/-l/ 

and we have the exponents 6 = and Ag = 2. The logarithmic divergence of the 
single-time correlator [25 3 reflects itself here in the logarithmic form of the scaling 
function. 

Finally, for al < 2 the same computation as for al > 2 goes through. Now, the 
exponent b = | — 1 is negative which means that the two-time autocorrelator 
diverges, in agreement with the earlier results for the equal-time correlators in ID 



r-dependence for s, t — s 1 

We use the asymptotic expression, valid for u ^> 1 and r 2 /u fixed 

(2^)1 ^ (~ 2^) 
which yields for arbitrary dimension ol, when introduced into (|17|) 



dz I ri (u) ■ . . . • I u (u) ~ — — -exp [ -— ) (23) 



(47T) 



2 



4 



(24) 



2 '2t + sJ \2 '2t-s. 
The incomplete Gamma function r(«, x) is defined by [3*9*] 

T(k,x) := / dte"*^ 1 (25) 

J X 

and has the following asymptotic behaviour for large or small arguments 

r(«, x) x«- x e~\ r(«, x) 0< ^ r(«) - — . (26) 

In the limit where both s and i — s become large, we recover eq. (fT9"j) as it should 
be. Furthermore, we explicitly see that the dynamical exponent z = 2. 
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1 2 3 4 5 6 5 10 15 20 25 



t/s r 2 /s 

Figure 2. Scaling plots of (a) the autocorrelation function G(0;t,s) and (b) the 
space-dependent correlation function G(r;t, s) for the critical bosonic contact process 
in three dimensions with apo = 1. In (b), the value of y = t/s = 2 was used. 



For illustration, we have also evaluated the integral (jl7)l numerically (with apo = 1). 
In Figure^ we compare the numerical results, for several values of s in three dimensions, 
with the analytical result eq. (fTT^l . We see that already for quite small values of s one 
has a nice data collapse which confirms the expected scaling behaviour. Furthermore, 
the agreement with the analytically calculated scaling function is perfect. In figure Eb 
we display the dependence on r, evaluated along the line r = (r, 0, . . . , 0). Again, the 
expected scaling behaviour is also confirmed and the curves agree with the analytical 
expression eq. (j2U). 

In the non-critical case, we have for the density, after rescaling t — > t/(2D) 

d . . 1 . . . 1 . . . _ uk — \t , 

— p(x, t) = -A x p{x, t) + -vp{x, t) with rj := — - — (27) 

This is easily solved and yields 

p(x,t) = poefr (28) 

if we choose again a homogenous initial distribution with mean density po- Depending 
on whether particle creation or annihilation dominates the density increases or decreases 
exponentially. Next, for the single-time correlator F(r,t) we use j2Hl eqs. (7,8)] which 
can be written after rescaling as (recall I = 1) 

^F(r, t) = AF(r, t) + V F(r, t) + a5 rfi pit) (29) 
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which is easily solved by introducing the particle-density (J28|) and performing a Fourier 
transform. This in turn allows to solve the equation of motion (JT4"Jl for the two-time 
correlator and we find 



F(r; t, s) = ple^ {t+s) + ap e^ {t+s) J dr e'^b (r; ~(t + s) - rj . 



(30) 



We consider the case r = and t and s in the ageing regime. As before, we use the 
asymptotic expression for 6(0, \ {t + s) — r) and find for the connected autocorrelator 

G(0;t,s) :=F(0;t,s)-p 2 e^ t+ ^ 



ap e^ t+s) 



.(31) 



( 47r )d/2 

Using the asymptotic behaviour eq. (J2T?|) for the Gamma-function for large arguments 
we obtain 

G(0; t, s) = -j^^ [(t - S )- d/2 exp - (t + s)^ exp + ,))]. (32) 

If rj is positive, then particle-creation outweighs particle-annihilation. The second 
term dominates and leads to an exponential divergence. On the other hand, if rj is 
negative, the first term involving e vt ^ 2 is the dominant one. At first sight, these results 
appear curious, since the leading exponential behaviour merely depends on t + s and t, 
respectively, and not on t — s, as might have been anticipated. 

A similar result had already been found in the inactive phase of the ordinary contact 
process |231 121] and we can understand the present result along similar lines. Consider 
the limits where \rj\ — > oo, such that diffusion plays virtually no role in comparison 
with the creation or annihilation processes. Then merely the creation and annihilation 
processes on a single site need to be considered. Correlators are given in terms of 
conditional probabilities and we now consider the two cases rj > and rj < 0. First, 
for 7] < 0, annihilation dominates and at late times there are only few particles left in 
the system. This is the same situation as in the inactive phase of the ordinary contact 
process. Then G(0; t, s) can only be non-vanishing if at time s a particle was present and 
it should only depend on t. On the other hand, for r\ > 0, the particle-density diverges 
exponentially and the number of possible reactions is conditioned by the density at time 
s, proportional to e vs ^ 2 , hence the dependence on t + s. Finally, the power-law prefactors 
relate to the diffusion between different sites. 



4. The bosonic critical pair-contact process 

For the bosonic pair-contact process, we have p = m = 2. The system (|T4^) of differential 
equations closes only for the critical case, i.e. for \i = p,k, and we shall restrict to this 
situation throughout. At criticality, the values of i and k do not influence the scaling 
behaviour. It was shown in [26J that in dimensions d > 2 there is a phase transition along 
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the critical line and we must therefore distinguish three cases, according to whether the 
reduced control parameter 

a' := (33) 
etc 

is negative, zero, or positive and where a was defined in (J 12)) and etc in (|13p. For d < 2 
one is always in the situation a' > 0. We recall the known results for the single-time 
autocorrelator F(0,t) which for large times behaves as 

• a < etc '■ 

F(0,t)%°° -4- (34) 



a = etc 



t—*oo 



^4$ ti- 1 for 2<d<4 



F(o,t)°£"t m-dm\a G (35) 



Po 



-t for d > 4 



4A2«c 

where A 2 is a known constant which is defined in |2fi| . 
a? > ftp or (i < 2 : 

F(0,t) Ap 2 exp{t/T ts ). (36) 

The known prefactor A and the time-scale r ts are dimension- dependent and positive. 
The exact expressions for them are not essential for our considerations and can be 
found in 



The solution of the equations of motion is quite analogous to the one of the bosonic 
contact process and the results from section 3 can be largely taken over. We find, again 
for initially uncorrelated particles of mean density po, 



F(r; t,s)=pl + a^dr F(0, r)b (r, ^(t + 



s)-r) (37) 



For t = s this formula agrees with eq. (21)] as it should. We are interested in 
the behaviour of the connected correlation function, see ()17)) . in the ageing regime. 
The analysis of eq. (|37jl is greatly simplified by recognizing that, quite in analogy with 
ageing in simple ferromagnets, there is some intermediate time-scale t p such that for 
times t < t p , one still is in some quasi-stationary regime while for t >t p one goes over 
into the ageing regime and that futhermore, the cross-over between these regimes occurs 
very rapidly jlSj- We denote by F a9e (0,r) the asymptotic ageing form of F(0,r) and 
write 

^drF(0,T)b(0,^(t + s)-r) 

= f "dr F(0, r)6(0, ht + s)-r)+s [ dv F age (0, sv)b(0, ht + s) - rv). (38) 

</0 1 Jt p /s 1 
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We denote the first term of the last line by Ci(t, s, t p ). Since we expect that t p ~ with 
< £ < 1 43,, we can replace the lower integration limit by in the second integral. 
This leaves us with the result 

G(0;t,s) = C 1 (t,s,t p ) + ^dTF age (0,r)b(o,^(t + s)-r^ . (39) 
On the other hand, we have the following rough estimate 



\Ci(t, s, t v )\ < t v max 

re[0,t p ] 



S2>1 i „/ vi d ( I 1 , . . t r 

« t p max F 0,r )\ s s 4vr -(t/s + 1) - 

T€[0,tp] \ \2 s 

In the three cases (i) a < ac and d > 2, (ii) a = «c and 2 < d < 4 and (iii) a = ac 
and d > 4 this leads by eqs. (|34|35)L respectively, to the upper bounds |Ci| < s^ d ^ 2 , 
s (C-i)rf/2 an( j s K-d/2 which vanish for s large more rapidly than G(0;t, s) ~ s 1_d / 2 , s° 
and s 2 ~ d / 2 , respectively and which are derived below. Hence Ci(t,s,t p ) is irrelevant for 
the determination of b and the scaling functions and will be dropped in what follows. 
Similarly, because of Ci(t,s,t p ) is non-leading if a > ac- 

We have also checked that for d > 4 this same result can be derived more explicitly 
using a Laplace transformation, along the lines of For the sake of brevity, these 
relatively straightforward calculations will not be reproduced here |41j . 



4-1. Ageing regime: r = and s, t — s ^> 1 

The most interesting cases are d > 2 and a < ac, which we will treat first. The 
asymptotic expression for F(0,t) is of the form F age (0,t) = Apffi, where £ and the 
prefactor A can be read of from equations (p?£ ]) -(|35 )) . 

We therefore get for the connected autocorrelator 



G(0;*, a ) = ^§4 fdrr* (\{t + s)-r 
(47r)2 Jo V 2 



d 



ap ° A ^-Jfito + ijVVj^t + iit + ai-^rl («) 



(f + l)(45r)3 \2 7 V2 y + 

where y = t/s and 2-^1 is a hypergeometric function. We deduce the general form of the 
scaling function 

fa(y) = - ^ ~ d (kv + 1)) ,F X (~ £ + 1; C + 2; . (42) 

(£ + i)(4tt)2 V 2 / V 2 y + iy 

and the exponents 

6 = -£ - 1 + ^ and A G = d (43) 
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and furthermore z = 2, see and below. For the different cases we obtain the following 
explicit expressions: 

• a < etc and d > 2: Here we have £ = and the pref actor is A = Therefore, 

we have a value of 

, d , x 

6=2-1- ( 44 ) 

The 2-Fi-function can be rewritten with the help of the relation (9.121,5) from 
so that we obtain an elementary expression for the scaling function: 



,,2 



fa(y) = - oN Uv + - (y - i)- f+1 ) • (45) 

a'(2ir) 2 [d — 2) v 7 

a = etc and d > 2: Here we have £ = | — 1 and £=lfor2<<i<4 and c? > 4 
respectively. This implies for 6 



, for 2<d<4 
b= < 1-2 for d>4 ' (46) 

d 

For 2 < d < 4 the scaling function is (with the prefactor A = — |2*o] ) 

1 1 \ 1 2" / I 

2l+Vn , v i fd d d 2 \ 



For ci > 4, the scaling function (|42|) can again be written as an elementary function 
with the help of a GauB recursion relation (eq. (9.137,4) from |4*2*]): 

faiv) ~ Pl 



4A 2 (27r)2(rf-2)(rf-4) 
x ((y + l)-f+ 2 - ( y - l)-f+ 2 + (d - 4)(y - 1)^ +1 ) . (48) 

a > olc or d < 2: Due to the exponential behaviour of F(0, r) we do not have a 
scaling behaviour in these cases. The integrals which enter the calculation are 
similar to those encountered in (|3*TJl and where the time-scale r ts and the factor A 
are defined in eq. (J3HJ): 

ap 2 Ae t+s /^ * 
(4yr) d / 



G(0; t, s) — — //f _ Nrf / 2 r ts 



d + i,*±i > )-rf-- + i, t " s 



(49) 



2 2r is y \ 2 2rt s 

Using the asymptotic behaviour of the the Gamma function ([26)1 . we see that the 
leading term in the scaling limit is 

G(0;t,s)^^(t- S y^e W ^. (50) 
(27r)2 r ts 

In contrast with the other cases treated before, the connected autocorrelator 
increases exponentially with the waiting time s. 
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4-2. r-dependence for s, t ^> 1 

In order to compute the r-dependence of the correlator, we follow the same strategy 
as in the last section. We use the approximation (|2*3|) which can be justified by an 
argument relying on an inequality similar to ([40)1 . We obtain the following results. 

• a < ac and d > 2: As we have -F(0, r) « —p\jct the computation is the same as 
for the contact process, compare equation (J23J. The result is 

-«Po /^ 2 V (i ~ 1} r/cf i r 2 \ /<* 1 r 2 \1 

G(r; ^ TjQ (t) [ r (2 - 2^) - r (2 - 5^)J '< 51 > 



10 



o 

o 



rn 



10 



10" 



10" 



-1 



10" 



" 1 


I , | , | , | , t 


T 


• s=25 




° S=50 = 




* s=75 




: s=100 




I analytical e 




I [eq.(44)] \ 




V 1 




\ : 


1 (a) 


- 1 l 


1 , 1 , 1 , 1 , -t 



1 2 



3 

t/s 



4 5 6 



O 10 



• s=25 
° s=50 

- s=180 

- s=250 

- analytical 

[eq.(50)] 




Figure 3. Scaling plots for the case a 1 < of (a) the autocorrelation function G(0; t, s) 
and (b) the space-dependent correlation function G(r; t, s) for the bosonic pair-contact 
process in five dimensions and with apo = 1. In (b), the value of y — t/s = 2 was used. 



a = etc and d > 4: We find the following result 
G(r;t,s) 



x 



4A 2 (47r)2 
t + s 



r 



-(1-1) 



1 



'2(t + s ; 

«2 



d 

2 '2t+s 



r 



r 



1 



d 

- -2, — 

2 ' 2(* 



'2(t-s) 

*.2 



It is straightforward to check consistency with (|4Hj) for the case s and t — 
larger than r 2 by using the asymptotic form (|26|) of the Gamma function. 



• (52) 
s much 
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• a = ac and 2 < d < 4: 

G(r;t,s) = - - NI P °. f'dTT^-H-a + s) -rHexp f — ; ^ . 

(d/2)|r(f-l)|7 2 ; ; ^\ 2 (t + s-2r)J 

We develop now the exponential function. The integrals are similar to those already 
seen so that we merely state the result 



„ d d d 2 \ 

x 2 Fi -+n, -;- + l;— . (53) 

2 1 V2 '2' 2 H/s + lJ K J 

a > ac or d < 2: Here again we develop the exponential function and obtain as 
final result 

oo 



1 I 1 I 1 t 




1 2 3 4 5 6 5 10 15 20 25 

t/s r 2 /s 

Figure 4. Scaling plots for the case a' = of (a) the autocorrelation function G(0; t, s) 
and (b) the space-dependent correlation function G(r; t, s) for the bosonic pair-contact 
process in five dimensions and with apo = 1. In (b), the value of y = t/ 's = 2 was used. 

In view of the numerous approximations needed to derive these results, it is of interest to 
check them numerically. In figure 01 we compare the results of the numerical integration 
of (|37j) with the analytical predictions (J44|45|51|) which apply for a < ac and d — 5. 
The nice collapse of the data shows that the scaling regime is already reached for the 
relatively small values of s used. The perfect agreement of the data with the analytical 
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results confirms that dropping the term C\ in is justified (and suggests that C\ 
should be considerably smaller than the rough estimate ([40))) . Similarly, we compare 
data for a = etc in 5D with the predictions ()48|52j) in figure 0] and similarly in 3D in 
figure Again the agreement is perfect. 



10 



-l 



t/5 



^lO" 2 



10" 



(a) 



• s=25 
□ s=50 



* s=75 
o s=100 
- analytical 

[eq. (46)] 



12 3 4 
t/s 



n 1 r 



5 6 



10" 



10 



-5 




(b) 



i 1 i 1 i 1 

• s=25 
□ s=40 

* s=60 
« s=75 
- analytical 

\\ [eq.(52)] 

\ 




5 10 15 20 25 
r /s 



Figure 5. Scaling behaviour of the two-time correlator G for the case a' = and 
three dimensions. The value apo was set to unity. In (b), the value of y = t/s = 2 was 
used. The data collapse occurs for 6 = 0. 



5. Response functions 

The response function of the first moment to an external field h(y, s) is given by 

8(a(x, t)) 



R(x,y;t,s) 



Sh(y,s) 
5.1. The contact process 



(55) 

h=0 



We apply the definition (|55|) on both sides of the equation of motion (|27|) and find, 
exploiting spatial translation-invariance, with r = x — y 

J^(r; t, s) = ^AR(r; t, s) + ^r]R(r; t, s) + 8(t - s). (56) 
This is the defining equation of a diffusion-type Green's function with the solution 

R(r; t, s) = r e^ s) & (r, ~(t - s^j G(t - s). (57) 
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bosonic contact process 


bosonic pair-contact process 


(X < Ctc 


ct = etc 


a 


2 


1 _ i 

2 x 


2 x 


b 


d _ 1 
2 L 


1 _ i 

2 x 


if 2 < d < 4 
f - 2 if d > 4 




d 


d 


d 


A G 


d 


d 


d 


z 


2 


2 


2 



Table 1. Ageing exponents of the critical bosonic contact and pair-contact processes in 
the different regimes. The results for the bosonic contact process hold for an arbitrary 
dimension d, but for the bosonic pair-contact process they only apply if d > 2, since 
a c = for d < 2. 



where b(r,t) was given in eq. (|T8j) and ro is a normalization constant. This expression 
is invariant under time-translations and does remain so even at criticality, see below. | 

5.2. The critical pair-contact process 

The equation of motion for the particle-density on the critical line does not change in 
comparison with the contact process, so that we can take over the result (JSTj) with r\ set 
to zero and have 

R(r; t, s) = r b (r, ~(t - s)J 6(t - s). (58) 

The autoresponse function in the scaling regime is obtained by setting r = and using 
the known asymptotic behaviour of the Bessel function, with the result (t > s) 

R{t, s) := R(0; t, s) ~ r (2n(t - s)y d/2 (59) 

from which we can read off the scaling function and the exponents a and \r 

«=|-1. fn(y) = j^(y-^ d/2 , A fl = d (60) 

We collect our results for the ageing exponents a, b, X G , \r, z in table [T] A few 
comments are now in order. First, for both the critical bosonic contact process and 
the critical bosonic pair-contact process with a < ac, we see by comparing the result 
for a with the corresponding ones for b, see table ^ that a = b. Together with the 

| Ageing is characterized by the existence of several competing stable stationary states (or a critical 
point) and time-translation invariance (tti) can no longer be requested. However, that does not mean 
that tti were always impossible and indeed tti can be recovered as a limit case, for certain specific 
values of the ageing exponents. A well-known example is the the response function of the spherical 
model quenched onto criticality (T = T c ) in d > 4 space dimensions 5 . 
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identity Xq = \r the critical ageing behaviour of these systems is quite analogous to the 
one of simple, reversible ferromagnets quenched to their critical temperature. Second, 
the critical bosonic pair-contact process with a = etc furnishes an analytically solved 
example where a and b are different. This is analogous to the result found for the ID 
and 2D critical ordinary contact process, where a = b — 1 was observed (2H1 IM| and 
where the relation A^ = Xr holds as well. However, there is no apparent simple and 
general relation between the exponents a and b for ageing systems without detailed 
balance. Third, our results for the critical bosonic pair-contact process provide further 
evidence against the generality of a recent proposal by Sastre et al. |44 to define a 
non-equilibrium temperature which was based on the implicit assumption that a = b 
would remain true even in the absence of detailed balance. Fourth, we can compare 
the form of the scaling function fii(y) of the autoresponse with the prediction of local 
scale- invariance quoted in eq. (J5|) . We find perfect agreement and identify a = a'. Fifth, 
we recall that for z = 2 there is a variant of local scale-invariance which takes the 
presence of a discrete lattice into account. It is possible to construct the corresponding 
representation of the Schrodinger Lie-algebra and then a response function transforming 
covariantly under it should read for t > s in al spatial dimensions 

i?(r;t, S )=r (t- S )^)/ 2 exp(^^)l r (^ , X r (u) := f[ I rj («) (61) 

where x is a scaling dimension and vq, M. are constants. Here the spatial distance r is an 
integer multiple of the lattice constant. Comparison with eqs. ()58I18|) shows complete 
agreement if we identify x — d/2 and M. = 1/2. 

6. Conclusions 

We have studied the ageing behaviour of the exactly solvable bosonic contact process 
and of the bosonic critical pair-contact process in order to get a better understanding 
on how the present scaling description of ageing, which is derived from the study of 
reversible systems with detailed balance, should be generalized for truly irreversible 
systems without detailed balance. This more general situation might be closer to what is 
going on in chemical or biological ageing than the reversible systems undergoing physical 
ageing, e.g. after a temperature quench. In comparison with the ordinary contact and 
pair-contact processes, these bosonic models permit an accumulation of many particles 
on a single site and this possibility does indeed affect the long-time behaviour of these 
models. Trivially, if either particle production or annihilation dominates, the mean 
occupation number will either diverge for large times or the population will die out, 
but if these rates are balanced there is a critical line where the mean particle-density is 
constant in time and the system's behaviour is more subtle. Indeed, on the critical line 
the long-time behaviour depends on how effectively single-particle diffusion is capable of 
homogenizing the system, see figure ^ For dimensions d < 2, there is always clustering 
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fn(y) 


f G (y) 


contact process 




{ y _ l)-f+l _ („ + 


pair 
contact 
process 


a < etc 


d>2 


(v-i)-* 


(„ _l)-f+i _ (y + 


a — ac 


2<d<4 


(v-i)-* 




d>4 


(y-i)-5 


(y + l)-§+ 2 - (y - l)-i+ 2 + (d - 4)(y - l)-§ +1 



Table 2. Scaling functions of the autoresponse and autocorrelation of the critical 
bosonic contact and bosonic pair-contact processes. They are only given up to a 
multiplicative factor, which may depend on the dimension. The logarithmic form ill'L'l) 
of fdy) f° r the 2D bosonic contact process may be obtained from a d — > 2 limit. 



at criticality, that is a few sites are highly populated and the others are empty. On the 
other hand, for d > 2 there is no clustering in the bosonic contact process, but in the 
bosonic pair-contact process there is a clustering transition at some a = ac such that 
clustering occurs for a > ac (where the diffusion is relatively weak) and there is a more 
or less homogeneous state for a < ac- 

This behaviour of the models also reflects itself in their ageing behaviour which we 
studied here. We anticipated in the ageing regime t, s 1 and t — s ^> 1 the scaling 
forms for the connected autocorrelator and autoresponse 

G(t, s) := G(0; t, s) = s- b f G (t/s) , R(t, s) := R(0; t, s) = s^f^t/a) (62) 

together with the asymptotics feniy) ~ y~ x a,R/ z as y ^> \ anc l our results for the 
exponents and the scaling functions are listed in tables ^ and EJ Specifically: 

(i) For d > 2, the ageing of the bosonic pair-contact process for a < ac lies in the same 
universality class as the bosonic contact process, since all critical exponents and 
the scaling functions coincide. Furthermore, the ageing behaviour in the bosonic 
contact and pair-contact processes does not depend on whether the parity of the 
total number of particles is conserved or not. All these systems have in common 
that their behaviour is strongly influenced by single-particle diffusion. One might 
wonder whether an analogy to the Janssen-Grassberger conjecture could be 
formulated. § 

(ii) While for d < 2, we still find a dynamical scaling behaviour in the critical bosonic 
contact process, there is no such scaling for the bosonic pair-contact process if 
a > ac, hence in particular for d < 2. Therefore, although both models have 
the same topology of their phase-diagrams for d < 2, see figure Hk, their ageing 
behaviour is different. 

§ An important ingredient of the models studied here seems to be that at criticality the mean particle- 
density stays constant. On the other hand, even if a 'soft' limit on the particle number per site is 
introduced, e.g. by a further reaction 3A — > 2A, the long-time behaviour is likely to be the one of the 
PCPD, as checked for the particle-density in |48| . 
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(iii) At the clustering transition a = etc in the critical bosonic pair-contact process, 
dynamical scaling occurs, but the ageing exponents a and b are different. Here 
the absence of detailed balance leads to a substantial modification of the scaling 
description with respect to what happens in critical ferromagnets. In particular, 
there is no non-trivial analogue of the limit fluctuation-dissipation ratio of critical 
ageing ferromagnets. A relation a / i, see (jfij). has also been observed in the 
ordinary critical contact process which also shares the property that = \r still 
holds |2l^ El] . However, according to the known examples, a simple and general 
relation between a and b does not seem to exist for systems without detailed balance. 
Further evidence from other non-equilbrium models would be welcome. 

(iv) On the other hand, the equality Xc = \r between the autocorrelation 
and autoresponse exponents, at the critical point of the steady-state and for 
uncorrelated initial states, seems to be a generic feature even for systems without 
detailed balance. 

(v) The form of the response function is in full agreement with local scale-invariance 
which confirms that the annihilation operator a(x) is a suitable candidate for a 
quasi-primary field || of local scale-invariance. We shall come back to a detailed 
analysis of the correlators from the point of view of local scale-invariance in a 
sequel paper. 

Explicit results were also derived for the space-dependent scaling functions of space- 
time correlator and responses. For the contact process, the space-dependent response 
function is given by eq. (JoTj) and the space-dependent correlation function by eq. 
For the critical pair-contact process, the space-dependent response function is given by 
eq. (|5H|). The space-dependent correlation function can be found in (i) eq. (|5Tj) for the 
case a < ac and d > 2, in (ii) eq. (J52|l for the case a = etc and d > 4 , in (iii) eq. 
for the case a = ac and 2 < d < 4 and in (iv) eq. (JHlj) for the case a > ac or d > 2 . 

Finally, we comment on a suggested relationship between the bosonic pair-contact 
process and the spherical model in the spherical model, a classical result by 

Berlin and Kac [23 states that the magnetization is spatially uniform, in particular the 
possibility that almost the entire macroscopic magnetization were carried by a single 
spin can be excluded. This is in remarkable contrast to the clustering transition which 
occurs in the bosonic pair-contact process. More formally, a closer inspection shows 
notable differences between the spherical constraint and the analogous equation used to 
derive the correlator F(0,t). This suggests that the analogies between the two models 
do not seem to have a deeper physical basis. 
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|| In conformal field-theory, a quasi-primary field transforms covariantly under the action of the 
conformal group |49|. This concept can be generalized to fields transforming covariantly under the 
action of a group of local scale-transformations, see JS] and references therein for details. 
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